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Abstract. Colliding winds of massive stars in binary systems are considered as candidate sites of high-energy non-thermal 
photon emission. They are already among the suggested counterparts for a few individual unidentified EGRET sources, but 
may constitute a detectable source population for the GLAST observatory. 

The present work investigates such population study of massive colliding wind systems at high-energy gamma-rays. Based 
on the recent detailed model (Reimer et al. 2006) for non-thermal photon production in prime candidate systems, we unveil the 
expected characteristics of this source class in the observables accessible at LAT energies. Combining the broadband emission 
model with the presently cataloged distribution of such systems and their individual parameters allows us to conclude on the 
expected maximum number of LAT-detections among massive stars in colliding wind binary systems. 
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INTRODUCTION 

Massive stars are hot, high-mass and luminous (~ 10 5 " 6 Lq) stellar systems that possess supersonic winds (V\y ~ 
several 1000 km/s) and accompagnied with a large mass loss rate (M ~ 10~ 6 "~ 4 MQ/yr). From more than a dozen 
galactic Wolf-Rayet (WR) binary systems non-thermal radio emission has been observed [e.g., 1], and this is in- 
terpreted as synchrotron radiation from a relativistic electron distribution. As a consequence inverse Compton (IC) 
scattering in the dense photospheric radiation fields of the stars as the presumably dominant leptonic emission pro- 
cess jjl foi is guaranteed to produce high-energy photons. Observational indications for massive stars systems being 
high-energy sources comes from population studies of the still unidentified EGRET sources: correlations with various 
massive star populations have been found [e.g.. l23ll . although their unambiguous identification must await instruments 
with significantly improved source localization capabilities (e.g. GLAST). The site of particle acceleration may be 
the instability driven shocks inside the winds B29I1 . shocked wind collision regions of multiple systems [8, 5] or the 
termination shock ll25ll . The following work concerns colliding winds (CWs) of WR-binary systems. Its main charac- 
teristics [8] has been meanwhile impressively confirmed by radio as well as X-ray observations [see e.g., for the 
nearby WR 147 system]: The collision of supersonic winds is creating a region of hot shocked gas that is separated by 
a contact discontinuity, and a forward and reverse shock follows. Its stagnation point is defined by the balance of the 
wind momenta (assuming homogeneous winds). The observation of synchrotron radiation indicates that the winds are 
permeated by magnetic fields supposedly originating from the massive star's surface. Estimates for surface magnetic 
field strengths range from below B s = 100 G (e.g. Oil) up to - 10 4 G [9]. In the following we fix this value to a reason- 
able 100 G, unless stated otherwise, and use the magnetic rotator theory (e.g. 12811 ') to estimate the field strength at the 
CWR. Typically >mG or higher field strengths are expected, assumed to be constant throughout the emission region. 
Note that because the stellar target photons for IC scattering arrive at the collision region from a preferred direction, 
the full angular dependence of the cross section has to be taken into account (e.g. Il22ll ). This leads to anisotropy effects 
like the emitted flux and cutoff energy being dependent on the sight line into the wind [see e.g.. l20ll . In binary systems 
the maximum flux level will be reached at phases where the WR- is behind the OB-star along the line-of-sight. Photon 
absorption due to yy pair production may affect the s pec trum at energies above ~ 70 GeV in CW systems. 

The most up-to-date catalogs of galactic WR- ll26ll27ll and OB-stars [ 10] contains more than 220 and 370 sources, 
respectively, with a detected binary fraction of ~ 40 — 50% Ir26ll . More than 30% are located in clusters and associ- 
ations. Its spatial distribution reflects the spiral arm structure of our Galaxy ([26] and references therein). The total 
number of WR-(systems) in the Galaxy has been estimated to ~ 8000 [17], total OB-star number counts in our Galaxy 
may reach values as high as 60000 [19]. This is at least an order of magnitude larger than the number of all presum- 



ably galactic, but still unidentified y-ray sources detected to date. How many of those massive WR-binary systems 
will GLAST's LAT be detecting at most? By combining a realistic physical source model with the statistics offered 
by the star catalogs we aim to answer this question in the following. The result will have implications for the chosen 
approach to probe the existence of WR-binaries in the the GLAST data. Whereas spatial coincidences are expected 
given the number of cataloged Wolf-Rayet stars and the anticipated number of y-ray source detections in the Galactic 
Plane above GLAST sensitivity, this study provides a prior for the size of a plausible gamma-ray bright WR-binary 
population as well as the individuals to be preferably investigated among the population. 

SAMPLE SELECTION AND PARAMETERS 

The construction of the sample for our study starts with the galactic WR-catalog ll26i l27ll . selecting all 88 binary 
systems. For more than a dozen of those sources the nominal shock location falls below the star's photosphere if wind- 
momentum balance determines the shock location, making the CW scheme an unsuitable setting for those. These are 
excluded in the following. Because of the distance squared dilution factor, close-by sources have a higher probability 
to be detected. Indeed, our statistical study shows that sources beyond ~ 3 kpc have a negligible chance to become a 
LAT source. Hence a distance cutoff in our sample was chosen at 4 kpc, which excludes a further 42 systems. Finally, 
a y-ray flux estimate on the basis of model simulations requires at least the knowledge of the shock location, estimated 
from its period or stellar separation and wind momenta of the participating stars, to provide reasonable estimates for 
the photon and field density at the shock. More than 10 sources have no period assigned, and hence are excluded from 
this study. We end up considering 21 galactic WR-binaries for potential LAT detectability. 

For estimating the maximum possible flux in LAT's energy range we consider in the following the IC component 
only as the most dominant y-ray production process, and use the Reimer et al. (2006) WR-binary model. The 
most important system parameters that impact the flux estimate significantly are the photospheric luminosity and 
temperature of the stars, their masses and momentum loss rates, their stellar separation and system distance. Only the 
distribution of the two latter ones shows an appreciable spread whereas all others have rather narrow ones, owing to 
the definition of a massive star. System parameters have been derived from j26l[Tll [Til |24ll . Because we're seeking the 
maximum number of likely LAT-sources among the WR-binary population, we allow the theoretically highest possible 
acceleration rate (E oc Ec/Rl(Vw/c) 2 , Rl = Larmor radius) and evaluate the y-ray output at a phase of maximum IC- 
flux. The system eccentricity e is known for only half a dozen sources in the considered sample with an average e « 0.2 
and the highest being e w 0.9. In the following we assume e = 0, and evaluate then the resulting y-ray flux uncertainty 
using extreme eccentricities. The same applies to the surface field, fixed at B s = 100 G now. An inclination angle of 90° 
is used for unknown system inclinations. Particle number and energy conservation arguments apply for the electron 
energy injection, as outlined in 112011 . An particle injection spectrum oc E~ 2 is used throughout this work. With these 
parameters we calculate the maximum flux expected at > 100 MeV for all sample sources using the model of I20I1 . 

RESULTS 

We consider a source to be detectable by the LAT if the maximum photon energy reaches at least the LAT energy 
range, and if the photon flux > lOOMeV is above LAT's sensitivity. For the latter we use an estimated factor 5 higher 
sensitivity level as published for LAT's high-latitude sensitivitjQ, 2 x 10~ 8 cm~ 2 s _1 , to account the fact that for the 
massive star systems being located preferentially in the Galactic Plane. Among the considered sample we find 6-7 WR- 
binaries to be detectable by the LAT in an one-year exposure, of which three are positional coincident with EGRET 
unidentified sources. All but one (for WR 70 only an upper limit in the radio band exists) turn out to be non-thermal 
radio emitters, and tend to be preferentially very-long-period binaries (see Fig. 1). The ratio of the wind momenta 
together with the stars' separation determines the distance of the shock jcob above the photosphere. The resulting 
distribution of the the shock-star separation (see Fig. 1) highlights the preference of LAT-detectible sources to possess 
larger values of jcob ■ This is plausible in the light of the severe IC-losses causing a low energy cutoff in the electron 
spectrum if the shock is located too deep in the dense photospheric radiation field, and thus inhibiting GeV photon 
production in shorter-period binaries. Finally, we find that only close-by systems (< lkpc) have a high chance to be 



http : //www-glast . s lac . Stanford . edu/sof tware/IS/glast_lat_perf ormance . htm 



12 3 4 

distonce [kpc] 

FIGURE 1. Normalized source number distribution of period (left), shock-star distance xob (middle) and source distance (right) 
of LAT-detectable WR-binaries (red/grey) from the source sample in comparison to all sample sources (red/grey and black). 

detected by the LAT, whereas WR-binaries located beyond ~ 1 kpc but up to ~ 3 kpc have a significantly reduced 
detection probability (see Fig. 1). 

To estimate the robustness of our result in view of unavoidable uncertainties in the parameter values, we vary each 
parameter (distance du stellar separation £>wr-ob, surface magnetic field B t , mass loss rates Mwr, wind velocities 
Voo,wr,ob> stellar temperature T e ff and luminosity Lbol) an d deduce the resulting number of LAT-detectable sources 
after each change. Specifically we alter period, di, B*, Mwr,ob> V°o,wr,oB' ^ eff anc ^ ^ bo1 a factor 2, 2, 10, 10, 2, 

0. 2 and 2 respectively. A variation of the system eccentricity to the extreme possible values adds a further factor 
±qi uncertainty to the stellar separation. As a result we find 6-7 systems in our sample are LAT-detectable with a 
total uncertainty of ±5, leading to the conclusion that the expected maximum number of detectable WR-binaries in a 
one-year GLAST-LAT exposure may be as high as eight, but could be as low as only one source from this population. 

Identification procedures may substancially benefit from physically motivated estimates of the expected source 
density of the respective population, like suggested here for the case of WR-binaries. Further parameter studies will 
include their spatial distribution, flux variability and orbital periodicity, spectral features, and signatures obtained at 
other wavelengths (see 1I2TI1 for a discussion applied to massive star systems). 
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